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Abstract. A subordinate Brownian motion A is a Levy process which can be obtained 
by replacing the time of the Brownian motion by an independent subordinator. In this 
paper, when the Laplace exponent (f> of the corresponding subordinator satisfies some 
mild conditions, we first prove the scale invariant boundary Harnack inequality for X on 
arbitrary open sets. Then we give an explicit form of sharp two-sided estimates on the 
Green functions of these subordinate Brownian motions in any bounded C 1 ' 1 open set. 
As a consequence, we prove the boundary Harnack inequality for X on any C 1 ' 1 open set 
with explicit decay rate. Unlike [KSV12b, KSV12c], our results cover geometric stable 
processes and relativistic geometric stable process, i.e. the cases when the subordinator 
has the Laplace exponent 

0(A) = log(l + A Q/2 ) (0<a<2,d>a) 

and 

0(A) = log(l + (A + m a/2 ) 2/a - m) (0 < a < 2, m > 0, d > 2) . 



1. Introduction 

Let d be a positive integer, let W = {Wt^x) be a Brownian motion in M. d starting at x 
and let S = (St : t > 0) be a subordinator independent of W, i.e. a Levy process taking 
values in [0, oo) and starting at 0. 

The Laplace exponent of a subordinator is a Bernstein function and hence has the repre- 
sentation 

<MA) = 6A+ f (1 - e~ xt ) n(dt) , (1.1) 

(0,oo) 

where b > and fj, is a measure on (0, oo) satisfying J (1 At) fi(dt) < oo, usually called 

(0,oo) 

the Levy measure of (p. If the measure [i has a completely monotone density, the Laplace 
exponent (ft is called a complete Bernstein function. 



2000 Mathematics Subject Classification. Primary 60J45, Secondary 60J75, 60G51. 

Key words and phrases, geometric stable process, Green function, Harnack inequality, Poisson kernel, 
harmonic function, potential, subordinator, subordinate Brownian motion. 

Research supported by Basic Science Research Program through the National Research Foundation of 
Korea(NRF) funded by the Ministry of Education, Science and Technology (0409-201 10087). 

Research supported in part by German Science Foundation DFG via IGK "Stochastics and real world 
models" and SFB 701. 

1 



2 



PANKI KIM AND ANTE MIMICA 



We define the subordinate Brownian motion X = (Xt,P x ) by Xt = Ws t - 

The aim of this paper is to obtain the following two-sided estimates of the Green function 
Gd{x, y) of X in a bounded C 1 ' 1 open set D C M. d in terms of the Laplace exponent <f> of 
the subordinator: 



\/0(fe(i)- 2 )*(<5D(y)- 2 ) j \* - i/l" +2 0(|z - sh 2 ) 2 ' 

where 5d(x) denotes the distance of the point x to D c and a A b := min{a, b}. Here and in 
the sequel, / x g means that the quotient ^ stays bounded between two positive numbers 
on their common domain of definition. 

The process X is, in particular, a rotationally symmetric Levy process. Recently there are 
many interests in studying potential theory of such processes. See, for example, [KMR, 
KSV12a, KSV12b, KSV12c, RSV06] and references therein. The purpose of this paper is 
to extend recent results in [KSV12b, KSV12c] by covering geometric stable processes and 
much more. 

Estimates of Green function for discontinuous Markov processes were first studied for ro- 
tationally symmetric a-stable processes in [CS98] and in [Kul97] independently. These 
results were extended later to relativistic a-stable processes and to sums of two inde- 
pendent stable processes in [Ryz02] and [CKS12] respectively. Recently, the first named 
author with R. Song and Z. Vondracek succeeded to obtain such estimates for a large class 
of subordinate Brownian motions in [KSV12b]. 

Still, the class considered in [KSV12b] does not include some interesting cases like geomet- 
ric stable processes or, more generally, the class of subordinate Brownian motions with 
Laplace exponent that varies slowly at infinity. Our approach covers a large class of such 
processes. 

Another feature of our approach is that it is unifying in the following sense: the sharp 
estimates of the Green function are given only in terms of the Laplace exponent (j) an d its 
derivative. 

Let us give a few examples of transient processes that are covered by our approach. 
Example 1 (Geometric stable processes) 

4>(\) = log(l + A^ 2 ), (0</?< 2, d>/3). 

Example 2 (Iterated geometric stable processes) 

01 (A) = log(l + A^ 2 ) (0</3<2) 
4> n+ i =(j) X o(j) n n e N, 

with an additional condition d > 2 1 ~ n (3 n . 
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Example 3 (Relativistic geometric stable processes) 

4>(X) = log (l + (a + m^ 2 ) 2/13 -rn\ (m > 0, < /3 < 2, d > 2). 

In order to obtain the sharp Green function estimates we first obtain the uniform boundary 
Harnack principle, with constant not depending on the open set itself. Such uniform 
boundary Harnack principle was first proved in [BKK08] and very recently generalized 
to a larger class of rotationally symmetric Levy processes in [KSV12c]. We adapt the 
approach in the latter paper in order to cover the class of subordinate Brownian motions 
with slowly varying Laplace exponents. Unlike the approach in [KSV12c], instead of the 
use of the Harnack inequality we use estimates of the Green function of balls near boudary 
obtained in [KM]. 

Further, our uniform boudary Harnack principle can be used to prove sharp Green function 
estimates for bounded C > open sets by adapting the method in [CKS12]. Even though 
we follow the roadmap in [CKS12], we needed to make significant changes due to the fact 
that now we do not have necessarily regularly varying Laplace exponents. 

To overcome such difficulties we use new types of estimates (not only in terms of the 
Laplace exponent itself, but also in terms of its derivative) of the jumping kernel and the 
potential kernel of the subordinate Brownian motions, which were obtained for the first 
time in [KM]. This type of estimates is essential in our approach. 

Let us be more precise now. In this paper we will always assume the following three 
conditions on the Laplace exponent 4> of the subordinator S: 

(A-l) (p is a complete Bernstein function; 

(A-2) the Levy density fj, of (f) is infinite, i.e. [i(0, oo) = oo ; 

(A-3) there exist constants a > 0, Ao > and 5 G (0, 1] such that 

< ax~ 5 for all x > 1 and A > A . (1.2) 

(f/{X) 



In the cases d < 2 and < 5 < \ we will sometimes further assume the following technical 
conditions: 

( A-4) If d < 2, we assume that the constant 5 in (A-3) satisfies d + 25 — 2 > and that 
there are ctq > and 

S e(l-i(l + l)A(2S + ^)) (1.3) 

such that 

4>'{^) . -So 



<P>(\) 



>a x~ d ° for all x > 1 and A > A . (1.4) 
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(A-5) If d > 2 and the constant S in (A-3) satisfies < S < ^, then we assume that there 
exist constants o\ > and Si G [5, 25) such that 

> ffi x 1 " 51 for all x > 1 and A > A . (1.5) 

0(A) 



Remark 1.1. (a) Note that (A-3) implies b = in (1.1), by letting x — > oo. 

(b) The condition (A-3) is implied by the following stronger condition 

Vx>0 lim = x^' 1 (0<a<2). (1.6) 

A^oo (j)'(X) 

In other words, (1.6) says that 0' varies regularly at infinity with index § — 1- A 
novelty here is the case a = 0. 

(c) The condition (A-4) is used only to obtain Green function estimates. 

Now we state the main result of this paper. By diam(D) we denote the diameter of D. 

Theorem 1.2. Suppose that X = (X t ,F x : t > 0,x G R d ) is a transient subordinate Brow- 
nian motion whose characteristic exponent is given by $(6) = <p{\6\ 2 ), 9 G R d satisfying 
(A-l)-(A-5). 

Then for every bounded C 1,1 open set D (see Definition 3.4) in M. d with characteristics 
(R,A), there exists c = c(diam(D),R,A,4>,d) > 1 such that the Green function Gr>(x,y) 
of X in D satisfies 

c~ l g D (x,y) < G D (x,y) < cg D (x,y) (1.7) 

with 

i i A, 0(1* -;/|- 2 ) <f(\z-y\-'') 

In particular, we extend the main result in [KSV12b]. 

Corollary 1.3. Suppose that X = (Xt : t > 0) is a subordinate Brownian motion whose 
characteristic exponent is given by 3>(0) = cf)(\9\ 2 ), 9 G M, d , where '• (0, oo) — > [0,oo) is a 
complete Bernstein function such that 

ci x a/2 < < c 2 x p/2 for all x>l and A > Ai . (1.9) 

0(A) 

for some constants ci, C2, Ai > 0, a, (3 G (0, 2) and a < (3. We further assume that 2/3 — a < 
1 ifd> 2 and /3 > 1, and that (A-5) hold with 5 = 1- /3/2. 

Then for every bounded C 1 ' 1 open set D in M. d with characteristics (R,A), there exists 
c = c(diam(D),R,A,<fi,d) > 1 such that the Green function Go(x,y) of X in D satisfies 
the following estimates: 

c^goix.y) < G D {x,y) < cg D (x,y) (1.10) 
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with 

- / \ /-, <fi(\ x -y\~ 2 ) 

g D (x,y) = 1A 



^</>(6 D (x)- 2 )<j>(8 D (y)- 2 ) I \x-y\ d <K\x-y\ 



-2) 



In [KSV12b], the above theorem is proved when, instead of (1.9), 4> satisfies 

0(A) x \ a / 2 £(\), A -> oo (0 < a < 2) (1.11) 
where £ varies slowly at infinity, i.e. 

Vx >0 lim^M = l. 

A^oo £{X) 

By Potter's theorem ([BGT87, Theorem 1.5.6(i)]), (1.11) clearly implies (1.9). 

Using Green function estimates, we prove the boundary Harnack principle for subordinate 
Brownian motions satisfying (A-l), (A-2), (A-3) and (A-5) in C 1, open set. We 
emphasize that in the next theorem we do not assume neither the transience nor (A- 
4). 

Theorem 1.4. Suppose that X = (X$,P Z : t > 0,x S M. d ) is a (not necessarily transient) 
subordinate Brownian motion satisfying (A-l), (A-2), (A-3) and (A-5) and that D is 
a (possibly unbounded) C 1 ' 1 open set in M. d with characteristics (R,A). Then there exists 
c = c(R, A,(f>) > such that for every r £ (0, ^p-], z € dD and any nonnegative function 
u in M. d that is harmonic in D n B{z,r) with respect to X and vanishes continuously on 
D c f]B(z,r), we have 



U ^^ C \^TfQ for every x,y^DnB { z^). (1.12) 



We remark that Theorem 1.4 cover the processes in Examples 1-3 without the assumptions 
on transience. 

Our paper is organized as follows. In Section 2 we record some preliminary results con- 
cerning subordinate Brownian motions obtained in [KM] . We start Section 3 by analyzing 
special harmonic functions in half-space and use these results to obtain key probabilistic 
estimates on C ' open sets. Section 4 contains estimates of Poisson kernel on balls which 
are used in Section 5 to obtain the uniform boundary Harnack principle on arbitrary open 
sets. After proving sharp Green function estimates in Lipschitz domains in Section 6, 
we finally obtain in Section 7 the boundary Harnack principle and sharp Green function 
estimates in C > open sets. 

Notation. Throughout the paper we use the notation f(r) x g(r), r — > a to denote that 
stays between two positive constants as r — > a. We say that / : R — > M is increasing if 
s < t implies f(s) < f(t) and analogously for a decreasing function. For a,i) 6 I, we set 
a A b := min{a, b}, a V b := max{a, 6}. For a Borel set A C we also use |A| to denote 
its Lebesgue measure. We will use ":=" to denote a definition, which is read as "is defined 
to be". 
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We will use the following conventions in this paper. The values of the constants Ci, C2, C3, 
C4 and £\ will remain the same throughout this paper, while c\, C2, . . . stand for constants 
whose values are unimportant and which may change from one appearance to another. 
All constants are positive finite numbers. The labeling of the constants c, c±, C2, . . . starts 
anew in the proof of each result. The dependence of the constant c on the dimension d 
will not be mentioned explicitly. 

2. Preliminaries 

By concavity, we see that every Bernstein function ift satisfies 

ip(tX) < Xift(t) for all A > 1, t > 0. (2.1) 
Thus A 1 — y ^j^- is decreasing, which implies 

A^'(A) < V(A) for all A > 0. (2.2) 

We first recall the following results from [KM]. 

Lemma 2.1. [KM, Lemma 3.1] Suppose that ift is a special Bernstein function, i.e., 
A 1 — V is also a Bernstein function. Then the functions rji,r)2: (0,oo) — > (0,oo) given 
by 

r ?1 (A) = AV(A) and % (A) = A 2 ||^ 
are increasing . □ 

Recall that we will always assume that the Laplace exponent (ft of S satisfies (A-l)-(A-3). 
We also recall the following elementary fact from [KM] which says that (A-3) controls 
the growth of (ft. 



Lemma 2.2. [KM, Lemma 2.2(h)] For every e > there exists c(e,a) > 1 such that 

< cx l ~ 5+£ for all x>\ and A > A ; (2.3) 



(ft(Xx) 



(ft(X) 

□ 



The analysis of 1-dimensional subordinate Brownian motions will be crucial in this ap- 
proach. Therefore we now consider an one-dimensional subordinate Brownian motion 
(Zt,P x ) with the characteristic exponent (ft{0 2 ). 

Let 

Z t ■= sup{0 V Z s : < s < t} 
be the supremum process of Z and let L = (Lt : t > 0) be a local time of Z — Z at 0. 
The right continuous inverse of L is a subordinator and it is called the ladder time 
process of Z. The process Z L -i is also a subordinator, called the ladder height process of 
Z. (For the basic properties of the ladder time and ladder height processes, we refer the 
reader to [Ber96, Chapter 6].) 
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Let k be the Laplace exponent of the ladder height process of Z. It follows from [Fri74, 
Corollary 9.7] that 

K(A) = exp jl/M« rfe j, VA>0. (2.4) 

By our assupmtions and [KSV12a, Propsoition 3.7] or [KMR, Proposition 2.1] we see that 
the ladder height process of Z has no drift and is not compound Poisson, and so the 
process Z does not creep upwards. Since Z is symmetric, we know that Z also does not 
creep downwards. 

Denote by V the potential measure of the ladder height process of Z. We will slightly abuse 
notation and use the same letter V to denote the renewal function of the ladder height 
process of Z, that is V(t) = V((0,t)). V is a smooth function by [KSV12a, Corollary 3.8]. 

Combining [KSV12a, Propsoition 3.7] and [Ber96, Proposition III. 1] the following result 
holds. 

Proposition 2.3. There exists a constant c > 1 such that for all r > 

; c ~ 1 < V(r) < . . 



□ 



We next consider multidimensional subordinate Brownian motions. Let W = (Wt 
(W/, . . . , Wf) : t > 0) be a Brownian motion in R d with 



E 



e ie-(w t -w ) 



e~ m \ V0eR<*,t>O, 



and let S be a subordinator independent of W with Laplace exponent <fi. In the remainder 
of this paper, we always assume that X = (Xt,W x ) is a subordinate process defined by 
Xt = Ws t - This process is a pure-jump symmetric Levy process with the characteristic 
exponent $(£) = <K|£| 2 ), i-e. 



E 

Moreover, $ has the representation 



HO = /(I - cos(C • y)Mdy) 



with the Levy measure of the form U(dx) = j(\x\) dx 

For any open set D, let us denote by td to denote the first exit time of D, i.e. 

t d = inf{« > : X f ^ D} . 

Using Proposition 2.3, the proof of the next result is the same as the one of [KSV12c, 
Proposition 3.2]. So we skip the proof. 
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Lemma 2.4. There exists c > such that for any r £ (0,oo) and xq £ M. d , 

^[^(,o,r)] < ^(OF(r - |x - xo|) X for x e B(x ,r). 



□ 



The process X has a transition density p(t, x, y) given by 

oo 

p(t,x,y) = J (4^)" d / 2 exp (-^)p(<? t G da). (2.5) 
o 

When X is transient, we can define the Green function (potential) by 

oo 

G(x,y) = g(\y - x\) = p(t,x,y)dt. 



o 

Note that g and j are decreasing. 
The following result is proved in [KM]. 

Proposition 2.5. Suppose (ft satisfies (A-l)-(A-4). Then we have 

j(r) X r _d_ y (r -2 ), r -> + . (2.6) 

If X is transient, then 

J^^r^Of. (2.7) 

□ 

As a consequence of (2.6) it follows that if <^> satisfies (A-l)-(A-4) then for any X > 0, 
there exists c = c(-fT) > 1 such that 

j(r)<cj(2r), Vr€(0,K). (2.8) 

The function j also enjoys the following property: if 4> satisfies (A-l)-(A-4) then there 
is a constant c > such that 

j(r + 1) < j(r) < cj(r + 1) for all r>l (2.9) 

(see [KM]). 

Let .D C be an open subset. The killed process X D is defined by 

X® = X t if t < t d and X® = A otherwise, 
where A is an extra point adjoined to D (usually called cemetery). 

The transition density of X D is given by 

p D (t,x,y) =p(t,x,y) -E x \p(t - t d , X TD ,y);T D < t] 

A subset D of M. d is said to be Greenian (for X) if X D is transient. When d > 3, any 
non-empty open set D C M. d is Greenian. An open set D C M. d is Greenian if and only if 



GREEN FUNCTION ESTIMATES FOR SBM 



9 



D c is non-polar for X (or equivalently, has positive capacity with respect to X). For any 

oo 

Greeninan open set D in M. d let Gd(x, y) = J PD(t, x, y) dt be the Green function of X D . 

o 

Gu{x,y) is symmetric and, for fixed y G D, Gi)(-,y) is harmonic (with respect to X) in 

D\{y}. 

The next two results are the key estimates in [KM]. 

Proposition 2.6. Suppose X is transient and 4> satisfies (A-l)-(A-4). There exist 
constants c\,C2 > and 61,02 G (0, ^), 2b\ < b 2 such that for all xq G K d and r G (0,1) 
we have 

C l 0(r-*) E i/ r B(xo,r-) < G B[x0jr) {x,y) < C 2 ^-2) %TB(so,r) 

/or a// x € B(xq, b\r) and y G B(xq, r) \ B(xq, b 2 r). □ 

Proposition 2.7. Suppose X is transient and <fi satisfies (A-l)-(A-4). There exist 
constants c\ > and a G (0, |) so £aa£ /or xo G K d and r G (0, 1) we aat;e 

E x[TB(x ,r)] > /or any x G -B(x , ar) . 

□ 

Before we state the Harnack inequality, we recall the definition of harmonic functions. 
Definition 2.8. Let D be an open subset o/M d . A function u defined on W 1 is said to be 

(i) harmonic in D with respect to X if 

E x [\u(X TB )\] < oo and u(x) = K x [u(X TB )] , x G B , 
for every open set B whose closure is a compact subset of D; 

(ii) regular harmonic in D with respect to X if it is harmonic in D with respect to X and 

u(x) = K x [u(X TD )] for any x G D . 

The following Harnack inequality is the main result of [KM]. 

Theorem 2.9 (Harnack inequality). There exists a constant c > such that for all 
xq G M rf and r G (0, 1) we have 

h{x\) < ch{x 2 ) for all x±, x 2 G B(xq, r/2) 

and for every non-negative function h: M. d — > [0, oo) which is harmonic in B(xo,r). □ 

By the result of Ikeda and Watanabe (see [IW62, Theorem 1]) the following formula is 
true 

W X (X TD EF) = J J G D (x,y)j(\z-y\)dydz (2.10) 

F D 
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for any F C D°. We define the Poisson kernel of the set D by 

K D (x, z) = J G D (x, y)j(\z - y\) dy, (2.11) 

D 

so that f x (X TD G F) = J K D (x, z) dz for any F C D°. 

F 

Proposition 2.10. Suppose X is transient and (ft satisfies (A-l)-(A-4). There exists 
c\ = ci(4>) > and C2 = C2(4>) > such that for every r G (0, 1) and xq G M. d , 

K B(xnr) (x,y) < cx-, ^y- XQ \-^ (2.12) 



y/<t>{r-l)<t>{{r-\x-xo\)- 2 ) 

^=35- 



< ^1%*^ (2.13) 



for all (x,y) £ B(xQ,r) x B(xo,r) and 



K BM (x ,y) > c 2 for all y G B(x ,r) C . (2.14) 

□ 



Proof. First using (2.8) and (2.9) to (2.11), then applying Lemma 2.4 and Proposition 
2.7, (2.12) and (2.14) follow easily (see the proof of [KSV12a, Proposition 4.10] for the 
details). (2.13) follows from (2.12) and the fact that 4> is increasing. □ 



3. Analysis on half-space and C 1,1 open sets 



In this section we establish key estimates which will be used in sections later in this paper. 

Recall that X = (Xt : t > 0) is the (i-dimensional subordinate Brownian motion defined 
by Xt = Ws t where W = (W 1 , . . . ,W d ) is a (not necessarily transient) <i-dimensional 
Brownian motion and S = (St : t > 0) an independent subordinator with the Laplace 
exponent (j) satisfying (A-l)-(A-3). In this section, we further assume that (A-4) holds. 

Let Z = (Zt : t > 0) be the one-dimensional subordinate Brownian motion defined by 
Zf.= W d St . 

Recall that V is denoted the renewal function of the ladder height process of Z. We use 
the notation 

K+ := {x = (xi, . . . , x d -i,x d ) := (x, x d ) G R d : x d > 0} 
for the half-space. 

Set w(x) := V((x d ) + ). Since Z t = Wg has a transition density, by using [Sil80, Theorem 
2], the proof of the next result is the same as the one of [KSV12b, Theorem 4.1]. We omit 
the proof. 

Theorem 3.1. The function w is harmonic in with respect to X and, for any r > 0, 
regular harmonic inR' 1 " 1 x (0,r) for X . □ 
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Using Theorem 3.1, (2.8) and (2.9), the proof of the next result is the same as the one of 
[KSV, Proposition 3.3]. 

Proposition 3.2. For all positive constants r$ and L, we have 

sup / w(y)j(\x -y\)dy < oo . 

x£R d :0<x d <L J 

B(x,r ) c nR^. 

□ 

Define an operator (A, 2) (A)) by 

Af(x) :=p.v. / (f(y)-f(x))j(\y-x\)dy 



lim / (f(y) - f(x))j(\y - x\)dy 

J :|x-j/|>e} 



J)(A):= i/:R d ^R:lim [ (f(y)-f(x))j{\y-x\)dy 

{ {i/6» d :|a:-i/|>e} 

exists and it is finite } . (3-1) 

Let Cq be the collection of C 2 functions in R d vanishing at infinity. It is well known that 
Cq C and that by the rotational symmetry of X, A restricted to Cq coincides with 

the infinitesimal generator C of the process X (see e.g. [Sat99, Theorem 31.5]). 

Since V is smooth by [KSV12a, Corollary 3.8], using our Theorem 3.1, (2.8) and (2.9), the 
proof of the next result is the same as [KSV, Proposition 3.3] or [KSV12b, Proposition 
4.2], so we skip the proof. 

Theorem 3.3. Aw(x) is well defined and Aw(x) = for all x £ W^_. □ 

In the rest of this section we aim to prove two key estimates of the exit probability and 
the exit time for C > open sets. Let us recall the definition of a C ' open set. 

Definition 3.4. An open set DinW d (d > 2) is said to be a C 1 ' 1 open set if there exist 
a localization radius R > and a constant A > such that for every z € dD, there exist 
a C 1 ' 1 -function ip = ip x : R ' -1 -)• R satisfying ip(Q) = 0, V^(0) = (0, . . . , 0), 

|| VVHoo < A, \Vi/)(x) -Vip(w)\ < A\x-w\, x.wef 1 " 1 

and an orthonormal coordinate system CS Z : y = (y\, ■ ■ ■ ,yd-i,yd) := (y 5 Vd) with origin 
at z such that 

B(z, R)C\D = {y = (y, y d ) £ B(0, R) in CS Z : y d > ip(y)}. 

The pair (R,A) is called the characteristics of the C 1 ' 1 open set D. By a C 1 ' 1 open set 
in M we mean an open set which can be expressed as the union of disjoint intervals so 
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that the minimum of the lengths of all these intervals is positive and the minimum of the 
distances between these intervals is positive. 

Remark 3.5. In some literature, the C 1, open set defined above is called a uniform C 1 ' 1 
open set since (R,A) is universal for all z £ dD. 

For x £ M. d , let Sgr>(x) denote the Euclidean distance between x and dD. Recall that for 
x £ ]R rf , 5d(x) is the Euclidean distance between x and dD. It is well known that any C 1 ' 1 
open set D with characteristics (R, A) there exists r\ > so that it satisfies: 

(i) uniform interior ball condition, i.e. for every x £ D with 8d{ x ) < r i there exists 
z x £ dD so that 

\x - z x \ = 5 dD (x) and B(x , n) C L>, 
for X0 = Zx + ril | E & T ; 

(ii) uniform exterior ball condition, i.e. for every y £ ~R d \ D with 8qd{u) < r i there 

5 aD (y) and % ,n)cK d \D, 



exists z y € so that 



for y = z + n 



y~ z v 

\v-*v\ ' 



Assume for the rest of this section that D is a C ' open set with characteristics (R, A) 
satisfying the uniform interior ball condition and the uniform exterior ball condition with 
the radius R < 1 (by choosing R smaller if necessary) . 



Lemma 3.6. Assume additionally that (A-5) holds. Fix Q € dD and let 

Kv) 



V(5 D (y)) y e B(Q, R) D D 
otherwise . 



There exists C\ = C\(A, R,4>) > independent of the point Q £ dD such that Ah is well 
defined in D n B(Q, ^) and 

\Ah(x)\<d for allxe DnB(Q,f). (3.2) 

Proof. We first note that when d = 1, the lemma follows from Proposition 3.2 and 
Theorem 3.3 by following the same proof as the one in [KSV12b, Lemma 4.4]. 

Assume now that d > 2. Fix x £ D(~)B(Q, -j) and let xq £ dD such that <5d(x) = jar — xo|- 

Denote by tp a C ' function and by CS = CS Xo an orthonormal coordinate system with 
xo chosen so that x = (0, x^) and 

B(x ,R)nD = {y = (y,y d ) in CS : y € B(0, R), y d > ^(y)} . 
We fix such ip and the coordinate system CS. 
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Define two auxiliary functions tpi,^' B(0, R) — > R by 

Mv) = R-VR 2 -\y\ 2 and = - (r - y/& 



- \y\ 2 



By the interior/exterior uniform ball conditions (with radius R) it follows that 

My) < i>(v) < My) for any yeDnB{x,%). (3.3) 

Now we define a function h x (y) = V(Sjj+(y)), where 

H+ = {y = (y,y d ) in CS: y d > 0} 
denote the half-space in CS. 

Since S H +(y) = (yd) + in CS, we can use Theorem 3.3 to deduce that 

Ah x (y)=0, Vy€H + . 

Now the idea is to show that A(h — h x ){x) is well defined and that there exists a constant 
Ci = Ci(A, R, 4>) > so that 

y |%)-MiOli(li/-*l)<fo<ci foran y e >°- ( 3 - 4 ) 

To do this we estimate the integral in (3.4) by the sum of the following three integrals: 
h= j (h(y)+h x (y))j(\y-x\)dy 

h = j (h(y) + h x (y))j(\y - x\) dy, where 



A 



A:={ye(DU H + ) n B{x, f ) : ^(y) < y d < ^(y)} 
h = J \Kv) ~ h x (y)\j(\y - x\) dy, where E := {y £ § ) : y d > ^i(y)} 

and prove that Ji + J 2 + I3 < C\ . 

To estimate l\ note that, by definition of h, h = on -B(Q> i?) c which gives 

h< sup / F(yd)i(|^-y|)dy + ci / j(\y\)dy < 00. 

< 2d < h e(^#) c £(o,#) c 

Here we have used Proposition 3.2 and the fact that the Levy measure is a finite measure 
away from the origin. 



11 
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Now we estimate 12- Denoting by rrid-i(dy) the surface measure, we obtain 



h < J J lA(y)(h x (y) + h(y))j [^/r 2 + \y d - x d \ 2 J m d M d y) dr 

\v\=r 



Since V is increasing and 



r- VR 2 -\y\ 2 < lJ i<\yl 

we can use (3.3) to deduce 

My) + %) < w (My) - My)) < w(2\y\) . 

Then, by the fact that j decreases, Proposition 2.3 and (2.6), we get 



_R 

4 



h < 2 



l A {y)V(2\y\)j(r)m d -i(dy)dr 



\y\=r 



R 

4 



< c 2 / r~ d - 2 7=4 = m d _i({y G ^4: |y| = r}) dr . 

VM ) 



Noting that (^(y) — (50 1 < = f° r |y| = r i we obtain 

"id-i({y: \y\ = r, My) < yd < My)}) < c 3^ d for r<f 

Thus, by the previous observation and the integration by parts we get 



R 
4 



J 2 < c 4 / r 



k 'Q=dr = cJr(- y ^))'dr 



VW^) 



< C4 



R 
4 



limr-^/ <^>(r 2 ) + / v </>(r 2 ) dr 





By Lemma 2.2 applied to a fixed e < 5 we see that there is a constant c = c(e) > so that 
which gives 



R 
4 



R 

4 



^2 < c 4 J ^(j){r- 2 )dr < c 5 /" r ~ 1+<5 " £ dr < oo . 
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In order to estimate I3, we consider two cases. First, if < y d = 5 R+ (y) < 5o(y), 

yd+R- x \y\ 2 
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h(y)-h x (y)<V(y d + R- 1 \y\')-V(y d )= J v(z)dz < R' L \y\'v{y d ), (3.5) 

Vd 

since v is decreasing. 

If y d = 5 H+ (y) > Sf)(y) and y E E, using the fact that 5d{v) is greater than or equal to 
the distance between y and the graph of ip\ and 



Vd 



R+^\y\ 2 + (R-y d ) 2 



< \y\ 2 < \vl 



y/\y\ 2 +{R-yd) 2 +{R-yd) ~ 2 ( R -Vd) ~ R ' 



we obtain 



Vd 



h x (y)-h(y)< J v(z)dz < R~ l \y\ 2 v(r~ ^\y\ 2 + (R- y d ) 2 ) . (3.6) 

R~^\y\ 2 +(R-yd) 2 

By (3.5) and (3.6), 

h<R^ J \y\ 2 v{Vd)j{\x - y\)dy 

En{y:y d <8 D (y)} 

+ R- 1 J \y\ 2 v (R - + (R- Vd) 2 ) 3(\x - y\)dy 

En{y.y d >8 D {y)} 
= :R-\L 1 +L 2 ). 



Since 



R 



Ec{z= (z,z d ) eR d : \z\ < j A ^2Rz d - zj and < z d < ^}, 



changing to polar coordinates for y and using(2.1), (2.2), (2.6) and Proposition 2.3, yields 



R 

2 



Li < c 6 J v(y d ) 




/ f ^2Ry d - y \ \ 

rU'((r 2 + \y d -x d \ 2 )- 1 ) 



(r 2 + \y d -x d \ 2 )( d + 2 )/ 2 



-dr 



dy d 



R 
2 



< c 7 



J v{y d ) 



\ ° 

/ T^V 2R yd-y 2 d \ 

r d 4>'{{r + \y d -x d \)- 2 )_ dr 



\ 



{r + \y d -x d \) d + 2 



dy d 
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If 5^ i by (A-3) 



/ 



R ^((r + \y d -x d \)- 2 ) dr 



(r + \yd ~ x d \) 2 
o 



2 f d>'((r + \y d -x d \)- 2 ) dr 

'J cp>((R + \y d - Xd \)-2)(r + \y d -x d \y 


1 \ -o 

2n 



+ 2 \ dr 



(i? + |y d -x d |)- 2 7 (r + |2/d-x d |) 2 
o 



=^((J?+|y d -x d |) 

<c 8 ^((ii+|y d -x d |)- 2 ) 

=c 8 0'((i? + |y d - s d |)- 2 )(12 + \y d - x d \r 25 J(r + \y d - x d \y 2+2S dr 

o 

<c^'((R + \y d - x d \Y 2 )(R + \y d - x d \y 2S \y d - x d \^ 28 ^ (3.7) 
Thus, in the case 6 > |, (3.7) implies 

2 

Li < c n J v(y d )dy d < c 12 V(§) < oo. (3.8) 
o 

If we apply Proposition 2.3 and use (A-5), we deduce 

Vd 

v(y d ) < - I v{s)ds = -V(y d ) < Cl3 l^((y d )- 2 )-V2 < 
Vd J Vd Vd 

o 

Therefore, in the case < S < 1/2 we use Si < 25 < 1, 5\ + 1 - 28 < 1 (8\ G [5, 2<5)) and 
the dominated convergence theorem to see that 

a 

2 

Xd ^ I WFW=^F^ dm (3 ' 9) 



is a strictly positive continuous function in Xd G [0, -j]. Hence it is bounded and, conse- 
quently, L\ < oo . 



If 5 = i, we can calculate 



Var+iw-^in^^iJ (310) 



o 



Since in this case 8i G [i, 1), the dominated convergence theorem 
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implies that 



x d H> 



(y d y 



111 



R 

\Vd-x d \ 



dyd 



(3.11) 



is a strictly positive continuous function in x d 6 [0, -|] and hence it is bounded. (3.10) 
and this imply that L\ < oo . 

Let us estimate L2. Switching to polar coordinates for y, and by the use of (2.6), Propo- 
sition 2.3 and (A-3) we get 



x d+ -j 



R ( y^F 7 



L2 < C15 



v(R - Vr 2 + (R-y d ) 2 )r d j((r 2 + \y d - x d \ 2 ) l ' 2 )dr 



dyd 



\ 

f / y/2Ry d -v* 



J 



< 



Cl6 



*d+% 



V 

/ V 2R ya-y 2 d 



v(R - ^r 2 + (R- y d r)4>'((r 2 + \y d - x^ 2 )- 1 ) 
(r 2 + \y d - x d | 2 )( rf + 2 )/ 2 



r dr 



dy d 



< C17 



v(R - ^r 2 + (R-y d )Z)<j>'((r + \y d - x d \)~ 2 ) ^ 



\ ° 



(r+ \yd ~ x d \) 2 



dy d . 



Since, for < r < R, 



and J2Ry d - yj < ^/R/2yJ2R - y d < R for < y d < x d + f , we have 



(yyft^Rvt-y* - r)) ^((r + |y d - x d |)~ 2 ) 



L2 " C17 / / (r + ^-„|) 2 ^ 
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Using (A-3), we see that with a := ^2Ry d — y 2 d and b := \y d — x d \, 
11 v(jy- d (a-r))^({r + b)- 2 ) 



I 



(r + bf 



-dr 



al s 2 v(y/y~d(a-r)W((r + br 2 ) ^ | j v(y/yj(a - r)) ^ 



(r + 6) 2 

o 



a/2 
Cl8 



(r + 6) 2 (6 + a/2) 2 (i^) 



(r + 6) 2 ( 1 - 5 ) 

v{y/y~d{a - r))dr 



a/2 



1-1 



(r + 6) 2 



a 



:= J B 1 (y d )+ J B 2 (y (i ). 



Using (A-5), Proposition 2.3, the estimate a > \Jy d R and the assumption 5\ — 5 < 5 



for < 6 < \, we deduce 



2 



B2(yd)dyd < c 20 < 



j\y d )- 5l+5 -^dy d <oo when0<5<± 
o 

R 

2 

V(R) J {y d ) 5 ~ 'Uy d < oo when \ < S < 1. 



By (3.7)-(3.11), we also have that 



I B 1 (y d )dy d < c 2 \ I v(y d ) / — - ; — ^ — drdy d < oo. 



o o 



+ \Vd - x d \) 2 



Thus L 2 < oo. 

Now we see that A(h — h x )(x) is well defined. Indeed, since h x (x) = h(x) and 

1 {yeDUH+: \y-x\>e}\Ky) ~ h x (y)\j(\y ~ x\) 

< 1 AuB( x ,^(Hy) + h x (y))j(\y - x\) + l E \h(y) - h x {y)\j{\y - x\) E L\R d ) 
we can use the dominated convergence theorem to deduce that limit 

lim / (h{y) - h x (y))j(\y -x\)dy 

elO J 

{yeDUH+: \y-x\>e} 

exists. Moreover, Ah{x) is then also well defined and satisfies |.A/i(x)| < C\ . 
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For a, b > 0, we define Dq(ci, b) := {y G D : a > Pq(u) > 0, \y\ < b}. 



Lemma 3.7. Assume additionally that (A-5) holds. There are constants R± = R±(R, A, <f>) G 
(0, ; R J a 

and x G Dq(t, r), 



(0, , R = ) and Ci = Ci(R,A,(f>) > 0, i = 1,2, such that for every r < R\, Q G dD 

16-y/ 1+(1+A) 2 



xy X TDQ(rr) ED) > Cl V(5 D (x)) (3.12) 

and 



E, 



r £>Q (r,r) 



< c 2 V(5 D (x)). (3.13) 



Proof. Without loss of generality we may assume that Q = and that tp-.W^ — > R is 
a C ' function such that in the coordinate system C5o 

5(0, R) H D = {{y, y d ) G J3(0, i?) in CS : y d > ■ 



The function p defined by p(y) = yd, — ip(y) satsifies 

p{y) 



< S D (y) < p(y) for all y G B(0, R) n D. (3.14) 



VI + A 2 
Define for a > 0, 

D tt = {?/GD:0< p(y) < a, \y\ < a} 

and a function 

V(£ D (y)) yeB(o,R)nD 



h(y) 



otherwise , 



Using Dynkin formula and he same approximation argument as in the proof of the Lemma 
4.5 in [KSV12b], from our Lemma 3.6 we have the following estimate for any open set 
U C B(0,f)nZ>: 

h(x) - C^tu < E x h{X TU ) < h{x) + C^tu (3.15) 
where C\ > is the constant from Lemma 3.6. 

By choosing A : = — , R we obtain 

J & Vl+(1+A) 2 

D r C Da C D(0, 1) D D for all r < ,4 . 
Indeed, y G Z? r and r > the following is true 

\y\ 2 = \y\ 2 + M 2 < r 2 + (\y d + IV#)I) 2 < (i + (i + A) V . (3.16) 

In particular, for r < A 

\y\ < ^1 + (1 + A)M = I . 

The idea is to choose Ao > 1 large enough so that (3.12) and (3.13) hold for r < Aq 1 A 
and x G D r . 
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We are going to show that there are constants c± , oi > such that for any A > 4 and 
x £ D^-iA the following two inequalities hold: 

E x [h(X TDx _ iA )} > d (V0(16A 2 ^ 2 ) - VHR- 2 )) ^xT Dx _ lA (3.17) 

X -n> x -i A G ^) > c 2 (0(16A 2 ^ 2 ) - 0(iT 2 )) ®xT Dx _ lA (3.18) 

Once we prove this, we can choose Ao > 4 so that 



0(16A 2 i?- 2 )>70(iF^+^. 

Then, for any A > Ao and x £ D^-ia we can use 

ci (V^IGA 2 ^- 2 ) - a/^TF 2 )) - Ci > d 

and (3.15) to get 

K(M*)) = /i(x) > E^X^^)] - CiE x t Da _ 1a > dE x r Dx _ lA , 

which proves (3.13) with Ri = X^A . 

Similarly, by (3.15), for any A > Ao and x £ D\-^a we have 
V(5 D (x)) = h{x) < K x [h(X TD ^_ iA )} + CxE x t Dx _ 1a 

< V(R)P X (x TD _ a eD)+ do, 1 (0(16A^~ 2 ) - ^(ir 2 ))" 1 !^.^, 



which together with (3.18) yields 

F(M*)) 



X Tn G D ) > 

X- 1 A 



V{R) + C4C2 1 (0(16A 2 fl- 2 ) - <A(i?- 2 )) 1 
This proves (3.12) with R\ = Ag l A . 

Now we prove (3.17). Note that for z G -Da -1 A anc ^ 2/ ^(0, 



-XR\ 
r > 



\z\ < X /1 + (1 + A 2 )A~ 1 A = A" 1 ! < 

implies 

j(\z-y\)>j(2\y\)>c 3 j(\y\)- 
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Then the Ikeda-Watanabe formula implies 

E x [h(X TDx _ lA )]> J J G D ^ lA (x,z)j(\z-y\)V(5 D (y))dzdy 

B(0,r)nD\D x _ lA D A _ lA 



>C 3 



B(0,R)nD\D 



G D ^ lA (x,z) dz 

B(0,R)nD\D x „ lA 

since < S D (y) by (3.14) . 

On the set E := {(y,y d ): 2A\y\ < y d , A _1 -f < \y\ < R} we have 



/ V(S D (y))j(\y\)dy 



\y\ < v / l + 4A 2 2/d and y d - j>(y) >y d - A\y\ > ^ . 
Since i£ C -B(0, J?) \ changing to polar coordinates gives 



R 



E x [h(X TDx _ iA )}> Ci E x [r Dx _ lA ] J j(r)F( i7T | 2F )r d - 1 dr 



with constant C4 > depending on A and d. 
Then (2.6) and Proposition 2.3 imply 

R 



r - 3 j^zX dr 



E x [h(X TDx _ iA )]>c 5 E x [r Dx _ lA ] J . v __ 5 _ 
We prove (3.18) similarly: 



X [X TD €D)>F X [X TD G B(0, J R)nD\B(0,A- 1 f 



R 



> cqE x [t Dx _ 1a ] J j{r)r d ~ l dr 



A-l R 



> c 7 E x [td x _ 1a ] J r- z 4>'{r- 2 )dr 

= 2- 1 c 7 E :c [t IVi J (0(16A 2 iT 2 ) - 0(i?" 2 )) . 
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□ 

4. Analysis of Poisson Kernel 

In this section we always assume that the Laplace exponent (f> of the subordinator S = 
(St - t> 0) satisfies (A-l)-(A-4) and corresponding subordinate Brownian motion X = 
(Xt,¥ x ) is transient. 

First we record an inequality. 

Lemma 4.1. For every Rq > 0, there exists a constant c(Ro,(j)) > such that 
A- 1 Ro 



X 2 / r 



1 4>'(r- 2 )dr+ j r- 3 (f>'(r- 2 )dr < c0(A 2 ), VA > ^. (4.1) 



A-i 



Proof. Assume A > ^. By (1.2), <p'(r~ 2 ) < Cl r 25 X 26 4>' (X 2 ) for r < A" 1 . Thus 
A-i Ro 
X 2 [ r~ l ct)'(r- 2 )dr+ [ r" 3 '•</>' '(r~ 2 )dr 



o a-i 



A-i i _ 2 Ro 

AV(A 2 ) / r - 1 ^—tdr-2 [ (<f>(r- 2 ))'dr 



0'(A 2 ) 

A-i 
A-i 

< c 2 ^(A 2 )A 2+25 J r- 1+2S dr + c 2 <t)(X 2 ) < c^' (X 2 )X 2 + <j>(X 2 )) < 2c 3 0(A 2 ) 


where we have used (2.2) in the last inequality. □ 
Recall that the infinitesimal generator C of X is given by 

Cf(x) = f (f(x + y)- f(x) - y ■ Vf(x)l { \ y \< e} ) j(\y\)dy (4.2) 

E'' 

for every e > and / G C 2 (M. d ) where C 2 is the collection of bounded C 2 functions in R rf . 
Using Lemma 4.1, we now prove [KSV12c, Lemma 4.2] under a weaker assumption. 

Lemma 4.2. There exists a constant c(<p) > such that for every f G C 2 (M. d ) with 
0</<l, 

\Cf r (x)\ < c<P(r~ 2 ) |l + sup J2\ 
w/iere / r (y) := f(^) and b := 2 f j(\z\)dz < oo. 

I*l>i 
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Proof. Set L x = sup ymd £ i>fc Then 

l/(« + y)-/W-yV/(«)| < 

which implies the following estimate 

\f r {z + y)- f r (z) - y ■ V fr(z)l{\ v \<r}\ < ^T^ 1 {|y|<'-} + 2 ' 1 {\y\>r} ■ 
Now, (2.6) and (4.1) yield 

\£fr(z)\ 

< J \fr(z + y) ~ fr{z) - y -V f r {z)l { \y\< r }\j(\y\)dy 

R d 

< IT / l{\ y \<r} lj £j(\y\)dy + 2 j l {r <\ y \< 1} j{\y\)dy + 2 f l {]v] > 1} j(\y\)dy 

R d R d R d 

< c^(r- 2 )(2 + ^) + 2 J j(\y\)dy, 

{M>i} 

where the constant c is independent of r £ (0, 1]. □ 

Lemma 4.3. For every a £ (0,1), there exists a positive constant c = c(a,<p) > such 
that for any r £ (0, 1) and any open set D with D C -6(0, r) 

¥ X (X TD £ -B(0,r) c ) < c(f>(r~ 2 )E x [T D ] for all x £ D n 5(0, or) . 

Proof. Using Lemma 4.2, the proof of the lemma is similar to that of [KSV12a, Lemma 
4.15]. We omit the details. □ 

Let A(x, a, b) := {y £ R d : a < \y — x\ < b} and recall that the Poisson kernel K B (x, z) of 
X in D is defined in (2.11). 

Unlike [KSV12c], instead of Harnack inequality we use Proposition 2.6 in the next propo- 
sition. 

Proposition 4.4. Let p £ (0,1). Then there exists a constant c((j),p) > such that for 
any r £ (0, 1) we have 

M 



1+P.r 
2 ' 



K B (o,s)(x,z)ds < cj^j{\z\) 

for all x £ -6(0, pr) and z £ A(0, r, r) . 

Proof. We split the Poisson kernel into two parts: 

K B (o,s)(x,z) = J G B{0:S) (x,y)j(\z - y\)dy = h(s) + I 2 (s) 

B(0,s) 
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where 

h(s)= J G B{0>s )(x,y)j(\z -y\)dy 

B(0,3s/4) 



h{s) = J G B{0iS) {x,y)j{\z - y\)dy. 
A(0,3s/V) 

First we consider I\(s). Since \z — y\ > -g\z\, we conclude from (2.7) that 

2s 

h(s)<j{jf) J G(x,y)dy< Cl j(\z\) J t^g^dt 



<c 2 



B(0,3s/4) 

mi 



Then, since \z\ < r, 

M \z\ 

ds 



h(s)ds<c 2 j(\z\) J 

1±E- i+P T . 

2 ' 2 ' 



On the other hand, by Proposition 2.6 and Lemma 2.4, 

-d-2<p'(s- 2 ) 



h(s)<c 3S - d ~'^l J Ey[T B{0 , s) }j(\z-y\)dy 

A(0,3s/4,s) 

< C4S ~d-2£^l f S(\z-V\) d 

A(0,3s/4,s) 
-d-2 <t>'{s-*) f j(\z-y\) 



- C4S ^- 2 ) 3/2 J y ffiZ3f*) dy ' 

A(0,3s/4,s) 

since s — \y\ < \z — y\. 

Observing that A(z, 3s/4, s) C -B(z, s) C A(0, \z\ — s, 2r) we arrive at 

-d-2 ^'(s- 2 ) f j(\v\) 



<K<>- 


2)3/2 


<j>'(. 


,-») 


4>(s- 


2)3/2 


<j>'(. 


,-») 


<t>(s~ 


2)3/2 



y v#H- a ) 

A(0,|z|-s,2r) 
2r 



\z\-s 



< ces- d - 2 J^VH(\z\-*))- 2 - 
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Then using the fact that s i-> 4>'{s~ 2 ) and s h-> (/)(s -2 ) -1 are increasing we obtain 

kl 1*1 

s 

</>(s-2)3/2 



I 2 (s)^<c 6 / ' m^wP V <K(M - a) 2 ) cfe 



2 ' 2 ' 

| z |-±±£ 

z l-2) 

0(|^| - 2 ) 3 / 2 



<c S tPr ]~^/l Zl ' 2) I s/W^dt. (4.3) 



o 



By Lemma 2.2 with e = | > for any a £ (0, 1) we have 

a a 


a 

< cW-^V'Ka" 2 ) / s~ 1+<5/2 ^s < c 8 ay/<t>(a- 2 ) 



o 



This together with (4.3) gives 

Izl / 1 _i_ \ - d ^ 2 



1+2 - 



(l+£ r ) d V(|z|-2) , 9 

< C 9 /(| 2 |-2)3/2 klV^CFl ) ^ C 10 j(\z\)^=rj. 



□ 



5. Uniform Boundary Harnack Principle 

In this section we give a proof of the uniform boundary Harnack principle for X in an 
arbitrary open set with the constant not depending on the open set itself. This type of 
the boundary Harnack principle was first obtained in [BKK08] for rotationally symmetric 
stable processes. Since, using results of previous section, the proofs in this section are 
almost identical to the one in [KSV12c, Section 5], we give details only on parts that 
require extra explanation. 

Recall that X = (X t , P x ) is a subordinate process defined by Xt — Ws t where W — 
(W*, P z ) is a Brownian motion in M. d independent of the subordinator S and the Laplace 
exponent (j) of the subordinator S satisfies (A-l)-(A-3). 

Using (2.8), (2.9), Proposition 2.10, Lemma 4.4 and the fact that for U C D 

K D (x,z)=Ku{x,z) + E x [K D {X Tu ,z)}, (x,z)eUxD c , (5.1) 
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the proof of the next result is the same as the one of [KSV12c, Lemma 5.2]. 

Lemma 5.1. Assume that X is transient and satisfies (A-l)-(A-4). For every p G (0, 1), 
there exists c = c((f>,p) > such that for every r G (0, 1), zq G M. d , U C B(zo,r) and for 
any (x,y) G (U H B(z ,pr)) x B(z ,r) c , 



( 



Ku(x,y) < c^-± 



\ 



j(\z - z \)Ku(z,y)dz + j(\y - z \) 



\u\b(z £ 



+p)r 



J 



The process X satisfies the hypothesis H in [SztOO]. Therefore, by [SztOO, Theorem 1], for 
a Lipschitz open set V C R d and an open subset U C V 



x (X Tu G dV) = and 



X {X TU G dz) = Ku(x,z)dz on V c . (5.2) 



Using (5.2) and Lemma 5.1, the proof of the next result is the same as the one of [KSV12c, 
Lemma 5.3]. 

Lemma 5.2. Assume that X is transient and satisfies (A-l)-(A-4). For every p G (0, 1), 
there exists c = c(<p,p) > such that for every r G (0, 1), for every zq G M. d , U C B(zo,r) 
and any nonnegative function u in M. d which is regular harmonic in U with respect to X 
and vanishes in U c n B(zo,r) we have 

u{x)<c-^T) j j(\y - z \)u(y)dy, x G U n B(z ,pr). 

(u\B(z ,£±fL))uB(z ,r)c 



We give a detailed proof of the next result. 

Lemma 5.3. Assume that X is transient and satisfies (A-l)-(A-4). There exists C2 = 
C2(d,4>) > 1 such that for every r G (0, 1), for every zq G M. d , U C B(zo,r) and for any 
(x,y) G UnB(z ,%)xB(z ,r) c , 



I 



VVB(*o,§) 



j(\z - z \)Ku(z,y)dz + j(\y - z \) 



< K u (x,y)<C 2 E x [T U ] 



j(\z - zo\)Ku(z,y)dz + j(\y - z \) . 



Vv*(*o,§) 



Proof. Without loss of generality, we assume z$ = 0. Fix r G (0, 1) and let U± := U D 
B(0, \r), U 2 := U n B(0, |r) and U 3 := U n B(0, |r). Let x G U n B(0, §), y G 5(0, rf. 
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By (5.1), 
Ku{x,y) 



E x [K u (X 7t , a ,y)] + K Ua (x,y) 
J K v (z, y)F x (X TU2 e dz) + K Ua (x, y) 



u\u 2 



J K u (z,y)F x (X 7V2 G dz) + / K v (z,y)K U2 {x, z)dz + K U2 (x,y) 
U 3 \U 2 U\U 3 

Ku(z,y)¥ x {X T edz)+ / K v {z,y) Gu 2 {x,w)j(\z - w\)dwdz 



u 3 \u 2 



u\u 3 



+ J Gu 2 (x,w)j(\y - w\)dw =: h + I 2 + h 



U 2 



From Lemma 4.3 and Lemma 5.1, we see that there exist c\ and c 2 such that 



h < ci ( sup Ku(z,y) ) 0(r )E x [r C / 2 ] < c 2 E x [t U2 



(5.3) 



Now using (2.8) and (2.9) one can check as in [KSV12c] that there exists C5 = cs(d, 0) > 1 
such that 

c^K x [t U2 ] f j(\z\)K u (z,y)dz<I 2 <c 5 K x [T U2 ] f j(\z\)K u (z,y)dz (5.4) 
u\u 3 u\u 3 



and 

c^E x [r U2 ]j(\y\) < h < c 5 E x [r U2 ]j(\y\) . 
The upper bound follows from (5.3)— (5.5). 
Using the strong Markov property, we get 



E x [ru] = E x [t U2 ]+E x ^ Xtu2 [tu] 



< E x [t U2 ] + ( supEM ) E X (X e B(0, f ) c ) 



(5.5) 



< E x [t U2 }+ c^{r- 2 )~ l <j){(^)- 2 )E x [T U2 } < c 7 E x [t U2 ], 

where in the second inequality we have used Lemma 2.4 and Lemma 4.3 and in last 
inequality we have used (2.1). 
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Since 



j(\z\)Ku(z,y)dz = J j(\z\)K u (z,y)dz+ J j(\z\)Ku(z,y)dz 
U\Ui U\U 3 U 3 \Ui 



< / j(\z\)Ku(z,y)dz + sup K v (z, y) 
J \zeu 3 
u\u 3 



j(\y\)dy, 



A(0,r/2,3r/4) 



by (2.6) and Lemma 5.1, 

j(\z\)Ku(z,y)dz 



U\Ui 

( 



< 



4 



1 + 



s~ 3 (j)'(s- 2 )ds 



\ 



j(\z\)K u (z,y)dz+j(\y\) 



\l\U; 



J 



4 



1 - 2- 



(r -2) 



{s- 2 ))'ds 



\ 



j(\z\)K u (z,y)dz+j(\y\) 



/ 



< 1 + C 9 



\ 



j(\z\)Ku(z,y)dz+j(\y\) 



Combining (2.1) and (5.4)— (5.6), we finish the proof of the lower bound. 



(5.6) 



□ 



Using Lemmas 5.2 and 5.3, the proof of the next result is the same as the one of [KSV12c, 
Lemma 5.5]. 

Lemma 5.4. Assume that X is transient and satisfies (A-l)-(A-4). For every zq £ M. d , 
every open set U C B(zo,r) and for any nonnegative function u in M. d which is regular 
harmonic in U with respect to X and vanishes a.e. on U c n B(zo,r) 

C^E x [tu] J j(\y-z \)u(y)dy<u(x)<C 2 E x [Tu} f j(\y - z Q \)u(y)dy 

for every x G U Pi B{zq, |) (where C2 is the constant from Lemma 5.3). 

As [KSV12c, Corollary 5.6], the last two lemmas immediately imply the following approx- 
imate factorization of the Poisson kernel. 

Corollary 5.5. Assume that X is transient and satisfies (A-l)-(A-4). Let zq £ M. d and 
D C M d be open. Then for every r G (0, 1) and all (x,y) G (DnB(z , §)) x (D c nB(z ,r) c ) 
it holds that 



C 2 1 E x [T DnB i yZOtr) \A D (zQ,r,y) < K D (x,y) < C 2 E x [T DnB ( Zo ^]A D (z , r, y) , (5.7) 



GREEN FUNCTION ESTIMATES FOR SBM 



29 



where 



A D (z ,r,y) := J j(\z - Zo\)K DnB ^ r )(z,y) dz 

(DnB(zo,r))\B(z ,%) 

+j(\y-*o\)+ J 3(\z-zo\)®z[K D {XT DnBlMOir) ,v) 



As the proof of [KSV12c, Theorem 1.1], Lemma 5.4 and (5.7) imply the uniform boundary 
Harnack principle with the constant not depending on the open set itself. That is why 
this type of result is called the uniform boundary Harnack principle. 

Theorem 5.6. Suppose that (j) satisfies (A-l)-(A-3). There exists a constant c = c((f>) > 
such that 

(i) For every zq G M. d , every open set D C M. d , every r G (0, 1) and for any nonnegative 
functions u,v in M. d which are regular harmonic in Dn B(zo,r) with respect to X 
and vanish a.e. on D c n B(zq, r), we have 

u(x) < u(y) 
v(x) ~ v{y) 

for all x,y G D Pi B(zo, |). 

(ii) If X is, additionally, transient and satisfies (A-4), then for every zq G M. d , every 
Greenian open set D C M. d , every r G (0, 1), we have 

KD(x 1 ,y 1 )K D (x 2 ,y2) < cK D (x 1 ,y 2 )K D (x 2 ,yi) 
for all x\,x 2 G D D B(zo, |) and all yi, y 2 G D° n B(zq, r) c . 

Proof. Under the assumption of transience and (A-l)-(A-4) the result follows from 
Lemma 5.4 and Corollary 5.5 (see proof of [KSV12c, Theorem 1.1]). 

If the process X is not transient, we can use argument similar as in the proof of [KM, 
Theorem 1.2, p. 17] where it is shown how to deduce Harnack inequality in dimensions 
d = 1,2 from Harnack inequality in dimension d > 3 (since in this case the process is 
always transient). Since we will use the argument in the proof of Theorem 1.4 again. Here 
we provide the detail for the readers' convenience. 

We use the notation x, = (x 1 , . . . , x d ~ l ) for x = (x 1 , . . . ,x d ~ 1 ,x d ) G M. d and X = 
^ S (x,x d ))- As in the proof of [KM, Theorem 1.2, p. 17], we have that for ev- 
ery x d G R, X = (X t , ¥%) is a (d — l)-dimensional subordinate Brownian motion with 
characteristic exponent $(£) = <^(|£| 2 ) f° r £ ^ 

Suppose (i) is true for for some d > 2 and let D be an open subset of M^" 1 and u, v : — > 
[0, oo) be functions that are regular harmonic in DC\B(xq, r) with respect to X and vanish 
on D c n B(xo,r) a.e. with respect to (d — l)-dimensional Lebesgue measure. 
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Let /: R d — > [0, oo) such that f(x,x d ) = u(x) and g: R d — > [0, oo) such that g(x,x d ) = 
v(x), which are regular harmonic in B(xq, r)xM with respect to X by the strong Markov 
property since T(B(x ,s)r\D)xR = hxf{i > : Xt ^ B(xq,s) n D}. Clearly they vanish on 
(B(xQ,r) x R) n (D x M) c a.e. with respect to d-dimensional Lebesgue measure. 

Thus, by applying the result to /, g, we see that there exists a constant c > such that 
for all x G open set D C and r E (0, 1) 



»ft)_/((5M))< c /p^ = c »M forall Sl326DnB(So , 



u(xi) g((xi,0)) " g((x 2 ,0)) v(x 2 ) 

Applying this argument first to d = 3 and then to d = 2 we finish the proof of the theorem. 

□ 



6. Green function estimates on bounded Lipschitz domain 

The purpose of this section is to establish sharp two-sided Green function estimates for X 
in any bounded Lipschitz domain D of M. d . 

Recall that we have assumed that X = (Xt,¥ x ) is the subordinate process defined by Xt = 
Wst where W = (W t , Pj;) is a Brownian motion in ~NL d independent of the subordinator S 
and the Laplace exponent <f> of the subordinator S satisfies (A-l)-(A-3). In this section 
we further assume that X is transient and that (A-4) also holds. 

We will first establish the interior estimates using Proposition 2.5 and Theorem 2.9. As 
in [KSV12b], once we have the interior estimates, we can apply Theorem 2.9 and the 
boundary Harnack principle (Theorem 5.6), and use the arguments of [BogOO, Han05] to 
get the full estimates for bounded Lipschitz domain D. 

Lemma 6.1. For every bounded domain D C M. d , there exists a constant C 2 = C 2 {d,(j), 
diam(D)) > such that 

G D (x,y) <c J X ~ V [. f 2 ^' { j X 2 ~ yl 2) forallx,yeD, (6.1) 
4>(\x-y\ l ) 1 

and for all x,y € D with b 2 1 \x — y\ < 5d(x) A 8d(v) 

r ( ^r-i \x-y\- d - 2 4>'{\x-y\- 2 ) , fi 9 , 

Gd(x,v) > C 3 flis,^ ( 6 - 2 ) 

where b 2 S (0, i) is the constant from Proposition 2.6. 

Proof. Since Grj(x,y) < g(\x — y\) and D is bounded, (6.1) is an immediate consequence 
of Proposition 2.5. 

Now we show (6.2). We have two cases: 
Case 1: \x — y\ < b 2 . 
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Since B(x, b% \x — y\)cD and y G A(x, \x — y\, b^ l \x — y\), we can use Proposition 2.6 to 
get 

G D (x,y) > G B(x b -i lx _ yl) {x,y) > c x ^b\\x- y \-*) E ^ T B(x,b^\x- y \)\ 



| X - i; |-d~2^ ( | x _ y |-2 ) 

<K\x-y\-*)* 



> c 2 — 



where in the last inequality we have used Proposition 2.7, (A-3) and the facts that b 2 £ 
(0, 5) and that the function r 1— » is decreasing. 

Case 2: \x — y\ > b 2 - 

In this case it follows that 5d(x) A Sd(v) > 1- Let xo £ 9B(y, 62)- Then 

62 1 |x — 2/| = 1 < <5d(x) A 5z)(y) 
and so, by the Case 1, we obtain 

G D (x ,y)>c 2 b \^-ip - (6.3) 



•'-1 



Since Gd(-, y) is harmonic in B(xq, y) U -B(x, y) (with respect to X), we can use Propo- 
sition 2.9 to deduce 

G D (x,y) = E x [G D (X TB ^ b2/4)1 y)]>E x [G D (X TB{ ^ 

> c 3 G D (x ,y)¥ x (X TB{x b2/4) e B(x Q , |)) . (6.4) 

By Proposition 2.10, (2.10) and (2.11) we get 

p *(^ s( ^ 2 /4) e s ^o, D) = y ^^o, *) 

>^y / -*!)**■ (6.5) 

Since \z — x\ < diam(L>), by the monotonicity of j we deduce 

to (v c d/ b 2 \\ ^ . bfj(diam(D)) 

Therefore, using (6.3)-(6.5) we deduce 

G D (x,y)>c 6 2 ^ 

C 6 b 2 \x-y\- d -'4/(.\x-y\- 3 ) 

- diam(D) <K\*-v\-' J F ' 

where in the last inequality we have used Lemma 2.1 (by considering the cases d = 1 and 
d > 2 separately) and the fact that 62 < \x — y\ < diam(l?) . □ 
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An open set D is said to be Lipschitz domain if there is a localization radius R\ > and 
a constant A > such that for every z G dD, there is a Lipschitz function (j) z : M^ -1 — > M. 
satisfying 

\4> z (x) - <t> z {w)\ < K\x - w\, 
and an orthonormal coordinate system CS Z with origin at z such that 

B(z,R 1 )nD = B(z,R 1 )n{y = (y,y d ) in CS Z : y d ><f> x {y)}. 
The pair (Ri, A) is called the characteristics of the Lipschitz domain D. 

Lemma 6.2. For every L > and bounded Lipschitz domain D with the characteristic 
A), domain D CW 1 there exists c = c(L, d, 4>, R\, A, diam(D)) > such that for every 
x,y G D with \x — y\ < L(5d(x) A Sniy)), 

G D (x, y) > c — -2^2 • (6.6) 

4>(\x-y\ l Y 

Proof. By symmetry of Gd we may assume Sd(x) < 5o{y)- Moreover, by Lemma 6.1 we 
can assume that L > 6 2 and so we only need to show (6.6) for 6 2 <5d(£) < |x — y| < L5d{x). 

Choose a point w G dB(x, &2^d(^))- Then Lemma 6.1 gives 

^^'^ " Cl ^feRFF ■ 

Since \y — w\ < \x — y\ + \x — w\ < (L + 1)5d(x) and Gd(x, • ) = Gd{ • ,%) is harmonic with 
respect to X in &2^d(^)) U B(w, &2£d0e))j using the assumption that D is a bounded 
Lipschitz domain, by Lemma 2.1 (by considering the cases d = 1 and d > 2 separately), 
Theorem 2.9 and Harnack chain argument obtain 

(M D (*))- d -V((& 2 Mx))- 2 ) 



G D (x,y) > c 2 G D (x,w) > c 3 - 



0((& 2 Mx)) 



-2^2 



|x — 2/| d 2 (//(|x - y| 2 ) 

0(k - y|~ 2 ) 2 

□ 

For the remainder of this section, we assume that D is a bounded Lipschitz domain with 
characteristics (i?i,A). 

Without loss of generality we may assume that R± < Since D is Lipschitz, there exists 
k = k(A) G (0, i) such that for each Q G dD and r G (0, R±), there exists a point 

A r (Q) G DnB(Q,r) satisfying fl(A r (Q),(cr) C DnB(Q,r). 



Recall that Gr>{-,y) is regular harmonic in D\B(y, e) for every e > and vanishes outside 
D. 

Fix zq € D with < 5d{zq) < R\ and set e\ := ■^ L . Define 

r(x,y) := 5 D (x)V 5 D (y)V \x - y\, x,y £ D 
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and 

B r x y \ . = / i A G D : 5 d( A ) > %r(x, y),\x-A\V\y-A\< 5r(x, y)} if r(x, y) < e x 

(6.7) 

Note that for every (x,y) € D x D with r(x,y) < E\ 

\8 D (A) < r(x,y) < 2k~ 1 5 d (A), AeB(x,y). (6.8) 

Set 

i/(( Sp(zo) \-2\ 
'Sp(zo)\-d-2 g U 2 ) > 



C A :=C 3 (1 Adiam(Z>))(! 



2 0((«sM)-2)2 



2 

By (6.1) and Lemma 2.1 (by considering the cases d = 1 and d >2 separately) we see that 

G D (x,z ) < C 4 for i;eD\%^). 

Now we define 

<? D (x) :=G d (x,z ) AC 4 . (6.9) 

We note that for 5d(z) < 6ei, 

g D {z) = G D (z,z ), 
since fei < an d thus |z - z \ > 6 D (z ) - 6ei > 

The following lemma follows from Theorem 2.9 and the standard Harnack chain argument: 

3 

Lemma 6.3. There exists c > 1 such that for every x G D satisfying 5d(x) > we 
have 

c _1 < go(x) < c. 

□ 

Theorem 6.4. Suppose X is transient and <p satisfies (A-l)-(A-4). If D is a bounded 
Lipschitz domain with characteristics (R\, A), then there exists c = c(diam(D), R\ , A, 0) > 
1 such that for every x,y £ D ; 

c?z)(^)^|x - y| d </H|x - y| 2 ) g D {Ay\x - y\ d (p{\x - y\ 2 ) 

(6.10) 

where go and B(x,y) are defined by (6.9) and (6.7) respectively. 

Proof. Since the proof is an adaptation of the proofs of [BogOO, Proposition 6] and 
[Han05, Theorem 2.4], we only give the proof when 6d(x) < 5d(jj) < § \x — y\. In this 
case, we have r(x,y) = \x — y\ 

By Theorem 2.9, we see that for all x,y G D and A\,A2 6 B(x, y), 

gr){Ai) is comparable to 5^(^.2) • 



> 

Q x ,Q y £dD with \Q X - x\ = 5 D (x) and \Q y - y\ = 5 D (y) . 



Set r = ' x 2 61 and choose 
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Pick points x\ = A Kr / 2 (Q x ) and yi = A Kr / 2 (Q y ) so that 

x, x\ E B(Q X , Kr/2) and y, y\ G B(Q y , nr/2) . 
Then one can easily check that \zq — Q x \ > nr and \y — Q x \ > r. 

Then Theorem 5.6 implies 

-i G D {xi,y) < G D (x,y) < G D {x u y) 
1 0d(zi) ~ 3£>(x) ~ £d(zi) 

for some ci > 1. 

Also, since |zo — Q y | > r and |xi — Q y \ > r, by Theorem 5.6 again, 
-l G D (x 1 ,yi) < G D (x 1: y) < G D (xi,yi) 
1 5r»(yi) ~ 9d{v) ~ 9d(vi) 

Therefore 

_ 2 . G D (x,y) 2 Gp(xi,yi) 
— r\ TT — c l" 



1 9D(xi)gD{yi) 9D(x)gD(y) (^1)3.0(2/1)' 

Now we can use Lemma 6.2 to get 

|si-yi|-"-y(|ii-yi|- 2 ) < G d( x ,V) < C2cj |gi-yi|- d -V(|xi-yi|- 2 ) 
gD{xi)g D {yi) </>(ki-yi|" 2 ) 2 - gp(x)gD(y) ~ fl£>(xi)si)(i/i) ^(ki-yil" 2 ) 2 

(6.11) 

for some c 2 > 1. 

Since < |xj — yi| < 2|x — y|, Lemma 2.1 (by considering the cases d = 1 and d > 2 
separately) yields 

|si-yi|- d -V(|gi-t/i|- 2 ) < 9 gd+2 |x-y|- d -V(9|g-y|- 2 ) < „ od+2 ^T^lVfe^Z!] 
^(Ixi-j/il- 2 ) 2 - ^ 0(9|x-3/|-^)^ ^ Z ' 6 <P(\x-y\-^ 



and 

|gi-»i|- d -V(|ai-yi|- 2 ) ■> o-l o-d-2 |x-y|- d - 2 ^(4- 1 |x-y|- 2 ) „_1 9 -d-2j£ z y|^V(i£-^f) 
<M|xi-?/l|- 2 ) 2 ~ 6 ' ^ 0(4-Mi-I/|- 2 ) 2 - 13 " Z 0(|a-y|- 2 ) a 

Therefore, 

2- d - 2 c- 1 c7 2 | a -y|- d - 2 / (|g-y|- 2 ) < G D (x,y) < 3 d + 2 c 2 c 2 ^-^f^-yT!) /g 

If r = then r(x, y) = |x — y| > e\ and so 

0i>(A) = 9d(z ) = C 4 and 5 D (xi) A 6 D (yi) > = ^ . 

Thus, in this case, Lemma 6.3 yields 

„-i ^ gp(A) 2 

c 3 < — p—r — < c 3 (6.13 
9D{xi)9D{yi) 

for some C3 > 1. 



GREEN FUNCTION ESTIMATES FOR SBM 



35 



In the case r < y we have r(x, y) = \x — y\ < e\ and r = \r(x, y). Hence 

S D (x 1 )AS D (y 1 )>^ = ^^. 

Since \xi — A\ V |yi — A\ < 5r(x, y) + \x\ — x\ + \y± — y\ < 5r(x, y) + 2nr < 6r(x, y), Theorem 
2.9 applied to gn gives 

-i s 9d(A) , _i 9d{A) 

c, < — -. — ■- < C4 and c 4 < — -. — r- < C4 (6.14) 
9d{xi) 9D{yi) 



for some constant C4 > 0. Combining (6.12)-(6.14), we get 

_-l 9D(x)gp(y) \x-y\- d - 2 <t>'(\x-y\- 2 ) < r , v < g D (x)g D (y) \ x -y\-d-2<f>i (\x-y\- 2 ) 

for all A e B(x,y). 



7. Explicit Green function estimates on bounded C ' -open sets 



The purpose of this section is to establish the explicit Green function estimates from 
Theorem 6.4 in the case of bounded C ' open sets. 

Theorem 7.1. Suppose that X = (Xt : t > 0) is a transient al- dimensional subordinate 
Brownian motion where the corresponding subordinator S has the Laplace exponent <f> 
satisfying (A-l)-(A-5). If D is a bounded C 1 ' 1 domain in M. d with C 1 ' 1 characteristics 
(R,A), then there exists c = c(R,A,(j), diam(D)) > such that 



(V{8 D (x)) Al) < g D (x) < c (V(5 D (x)) A 1) for all x G D. 



(7.1) 



Proof. The proof follows from the proof of [KSV12b, Theorem 4.6] by using our Propo- 
sition 2.5, Lemma 3.7 and Theorem 5.6 . □ 



Proof of Theorem 1.2. When D is connected, Theorem 1.2 follows from [KSV12b, 
(4.38)] and our Theorems 6.4 and 7.1. 

Next we assume that D is not connected. The proof below is similar to the one in 
[CKSV12]. 

Let (R, A) be the C > characteristics of D. Note that D has only finitely many components 
and the distance between any two distinct components of D is at least R > 0. 

Assume first that x and y are in two distinct components of D. Let D(x) be the component 
of D that contains x. Then by the strong Markov property and (2.10) we obtain 



Gd (x, y) = E x . G D {X TD(x) , y) 



T D(x) 



I 



j(\X s - z\)G D (z,y)dz ds 



\p\D(x) 
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Consequently, 

j(dmm(D))E x [r D{x) ] J G D (y, z)dz < G D (x,y) < j(R)E x [r D{x) ] f G D (y,z)dz. 

D\D(x) D\D(x) 

(7.2) 

Applying the two-sided estimates (1.7) established in the first part of this proof to D(x), 
we get 



By (7.3) we get 



< E x [t D (x)) < 



G D (y,z)dz > J G D(y) {y,z)dz = E y [r D[y) ] > ^ ( ^ 2 fa) _ 2) - 

D\D(x) D(y) 



(7.3) 



On the other hand, (2.10) and (7.3) imply 



G D (y,z)dz < E y [T D ] = E y T D ( y ) 



+ E. L 



E x \t d ] 



D\D{x) 



< 



< 



C3 



+ E, 



T D(y) 



j(\X s - z\)E z [r D ]dzds 

D\D(y) 
+ C i j{R)E y [T D ( y )} < 



C5 



V^Dfe)" 2 ) 



We conclude from the last three displays and (7.2) that there is a constant cq > 1 such 
that 



< G D (x,y) < 



(7.4) 



Noting that 

R < \x — y\ < diam(D) 
when x and y are in different components od D, we obtain (1.7). 

Now we assume that x,y are in the same component U of D. Applying (1.7) to U we get 



G D (x,y) > Gu(x,y) > c 7 



1A <P(\x-y\- 2 ) ) \x-y\- d -^'{\x-y\-l) 

V0(%W- 2 W5Lr(j/)- 2 ) J ^-y\-' 2 )' 2 



c 7 1 A 



<t>(\x-y\~ 2 ) 



\x-y\- d -^'(\x-y\- 2 ) 

— ^rra ■ 
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For the upper bound, we use the strong Markov property, (2.10) and (7.3)-(7.4) to get 

G D (x,y) 
= G u (x,y)+E x [G D (X Tu ,y)] 

TV 



< ro ( 1 A <K\*-V\~ 2 ) ) \x-v\- d - 2 #(\*-v\- 2 ) , W 



D\U 



j(\X s - z\)G D (z,y)dzds 

D\U 

(\ c 9 / — dz ^ z 

1A ^(|a-y|- 2 ) ] \x-y\- d - 2 <P'(\x-y\- 2 ) d\u , ?5 ^ 

V#Ms)- 2 )<K<W!/)- 2 ) / Ms-vl" 2 ) 2 \Z*OM*)- 2 )0(<My)- 2 ) ' 1 ' j 



< c 8 

Since Z? is bounded, we get 

i f dz < \D\ 

y/^(S D (x)-^(5 D {y)-^) J y/<K6 D (z)-2) - V^(^)- 2 )*(<5Dfe)- 2 )0(diam(D)-2) 
D\U 



V*(«5r>W- 2 )*(«5r>(!/)- 2 ) 



which together with (7.5) gives 



G D (x,y) < en 1 A 



^(l^-y|- 2 ) 



□ 



Before we prove Corollary 1.3, we record a simple fact. 

Lemma 7.2. 7/(5* G (0, 1) and is a Bernstein function satisfying 

> ex 1 " 5 * /or a// x > 1 and A > Aq , (7-6) 



V>(Ax) _ 



/or some o~ > 0. T/ien i/tere exists a constant c > suc/i i/iaf f/H^) < c\ip'{X) for all 
A > A . 

Proof. Let a\ = 2 V (-) . Since ■?// is decreasing, 

aiA 

(oi - l)A^'(A) > J if/(t)dt = tpiatX) - (7.7) 



A 



Since ^(aiA) > cra^ *V(A) by (7.6), from (7.7), we get 

(oi - l)A^'(A) > (aa^ 5 * - l)^(A) > ij>(\). 

□ 
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Proof of Corollary 1.3. It follows from the assumptions of the corollary that we can 
apply Lemma 7.2. Thus by (2.2) and Lemma 7.2 we obtain < </>(A) < c\(j)'{\) for 

all A > A . Therefore (A-l)-(A-5) hold and (1.7) is equivalent to (1.10). □ 

Proof of Theorem 1.4. When d = 1, the theorem follows from Proposition 2.3, Theorem 
3.1 and Theorem 5.6 (i). 

Note that the result in [CKSV12, Lemma 4.2] is true in our case too. By this result, 
Theorem 2.9, Theorem 5.6 (i) and Theorem 1.2, the proof of Theorem 1.4 is the same as 
the proof of [KSV12b, Theorem 1.3] when d > 3. 

Note that if D is a C 1 ' 1 open set in R rf_1 with characteristics (R, A), then Dxlis clearly 
C 1,1 open set in W d with the same characteristics (R, A). Thus the case d = 2 can be 
handled in the same way as in the proof of Theorem 5.6 (i) . □ 

References 

[Ber96] J. Bertoin, Levy processes, Cambridge University Press, Cambridge, 1996. 

[BGT87] N. H. Bingham, C. M. Goldie, and J. L. Teugels, Regular variation, Cambridge University Press, 
Cambridge, 1987. 

[BKK08] K. Bogdan, T. Kulczycki, and M. Kwasnicki, Estimates and structure of a-harmonic functions, 

Probab. Theory Related Fields 140 (2008), no. 3-4, 345-381. 
[BogOO] Krzysztof Bogdan, Sharp estimates for the Green function in Lipschitz domains, J. Math. Anal. 

Appl. 243 (2000), no. 2, 326-337. 
[CKS12] Z.-Q. Chen, P. Kim, and R. Song, Dirichlet heat kernel estimates for A Q//2 + A 13 / 2 , to appear 

in Illinois J. Math. (2012). 

[CKSV12] Z.-Q. Chen, P. Kim, R. Song, and Z. Vondracek, Boundary Harnack principle for A + A a//2 , 

Trans. Amer. Math. Soc (2012). 
[CS98] Z.-Q. Chen and R. Song, Estimates on Green functions and Poisson kernels for symmetric 

stable processes, Math. Ann. 312 (1998), no. 3, 465-501. MR 1654824 (2000b:60179) 
[Fri74] Bert Fristedt, Sample functions of stochastic processes with stationary, independent increments, 

Advances in probability and related topics, Vol. 3, Dekker, New York, 1974, pp. 241-396. 
[Han05] Wolfhard Hansen, Uniform boundary Harnack principle and generalized triangle property, J. 

Funct. Anal. 226 (2005), no. 2, 452-484. 
[IW62] N. Ikeda and S. Watanabe, On some relations between the harmonic measure and the Levy 

measure for a certain class of Markov processes, J. Math. Kyoto Univ. 2 (1962), 79-95. 
[KM] P. Kim and A. Mimica, Harnack inequalities for subordinate Brownian motions, (preprint, 

2011). 

[KMR] M. Kwasnicki, J. Malecki, and M. Ryznar, Suprema of Levy processes, Ann. Probab., (Preprint, 
2011). 

[KSV] P. Kim, R. Song, and Z. Vondracek, Potential theory of subordinate Brownian motions with 

Gaussian components, preprint(2012). 
[KSV12a] , Potential theory for subordinate Brownian motions revisited, Stochastic Analysis and its 

applications to Mathematical Finance, Interdisciplinary Mathematical Sciences, vol. 13, World 

Scientific, 2012. 

[KSV12b] , Two-sided Green function estimates for killed subordinate Brownian motions, Proc. 

London Math. Soc. 104 (2012), 927-958. 

[KSV12c] , Uniform boundary Harnack principle for rotationally symmetric Levy processes in gen- 
eral open sets, to appear in Science in China (2012). 

[Kul97] T. Kulczycki, Properties of Green function of symmetric stable processes, Probab. Math. Statist. 
17 (1997), no. 2, Acta Univ. Wratislav. No. 2029, 339-364. MR 1490808 (98m:60119) 



GREEN FUNCTION ESTIMATES FOR SBM 



39 



[RSV06] M. Rao, R. Song, and Z. Vondracek, Green function estimates and Harnack inequality for 

subordinate Brownian motions, Potential Anal. 25 (2006), 1-27. 
[Ryz02] M. Ryznar, Estimates of Green function for relativistic a-stable process, Potential Anal. 17 

(2002), no. 1, 1-23. MR 1906405 (2003f:60087) 
[Sat99] K.-I. Sato, Levy processes and infinitely divisible distributions, Cambridge University Press, 

Cambridge, 1999. 

[Sil80] Martin L. Silverstein, Classification of coharmonic and coinvariant functions for a Levy process, 

Ann. Probab. 8 (1980), no. 3, 539-575. 
[SztOO] P. Sztonyk, On harmonic measure for Levy processes, Probab. Math. Statist. 20 (2000), 383- 

390. 

Department of Mathematical Sciences and Research Institute of Mathematics, Seoul Na- 
tional University, Building 27, 1 Gwanak-ro, Gwanak-gu Seoul 151-747, Republic of Korea 

E-mail address: pkim@snu.ac.kr 

Fakultat fur Mathematik, Universitat Bielefeld, Postfach 100131, D-33501 Bielefeld, Ger- 
many 



E-mail address: amimica@math.uni-bielefeld.de 



